A perturbative approach to a class of Fokker-Planck equations 
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In this paper we present a direct perturbative method to solving certain Fokker-Planck equations, 
which have constant diffusion coefficients and some small parameters in the drift coefficients. The 
method makes use of the connection between the Fokker-Planck and Schrodinger equations. Two 
examples are used to illustrate the method. In the first example the drift coefficient depends only 
on time but not on space. In the second example we consider the Uhlenbeck-Ornstein process with 
a small drift coefficient. These examples show that the such perturbative approach can be a useful 
tool to obtain approximate solutions of Fokker-Planck equations with constant diffusion coefficients. 
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I. INTRODUCTION 



The Fokker-Planck (FP) equation is one of the basic tools used to deal with fluctuations in various kinds of systems 
^\ • It has found applications in such diverse areas as physics, astrophysics, chemistry, biology, finance, etc. Owing 

to its wide applicability, various methods of finding exact and approximate solutions of the FP equations have been 
developed. Most of the methods, however, are concerned with FP equations with time-independent diffusion and drift 
(— I ' coefficients [ij. Generally, it is not easy to find solutions of FP equations with time-dependent diffusion and drift 
rS |. coefficients. 

One of the methods of solving FP equation with time-independent diffusion and drift coefficients is to transform the 
FP equation into a Schrodinger-like equation, and then solve the eigenvalue problem of the latter. The transformation 
• to the Schrodinger equation of a FP equation eliminates the first order spatial derivative in the FP operator and creates 
M ' a Hermitian spatial differential operator. This method is useful when the associated Schrodinger equation is exactly 
solvable; for example with infinite square well, harmonic oscillator potentials, etc. Several FP equations have been 
exactly solved in this way [ij . The method can also be useful for approximate results if a finite part of the spectrum 
of the associated Schrodinger equation can be exactly known ^2*1. When this happens the Schrodinger equation, and 
QQ ' hence the associated FP equation, are called quasi-exactly solvable _3J. 

\ In this paper we would like to show that the connection between the FP and Schrodinger equations can also be 
' useful even when the drift coefficients are time-dependent. Based on this connection, we present a direct perturbative 
CN '\ method which can be used to solve a certain class of FP equations. 
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II. PERTURBATIVE APPROACH TO FP EQUATIONS 



. „ , In one dimension, the FP equation of the probability density P{x^ t) is 1] 

^ ' §-Wix,t)=^-^D^'\x,t) + ^D('\x,t)^W{x,t). (1) 

The functions D^-^''{x,t) and Z)(^^(x,t) in the FP equation are, respectively, the drift and the diffusion coefficient. 
The drift coefficient represents the external force acting on the particle, while the diffusion coefficient accounts for 
the effect of fluctuation. The drift coefficient is usually expressed in terms of a drift potential U{x,t) according to 
L»(i)(x,t) = -dU{x,t)/dx. In this paper we shall be concerned with an important class of FP equations, namely, 
those with constant diffusion coefficients D^^''{x,t) — D > 0. 

The FP equation is closely related to the Schrodinger equation. A FP equation with time-independent drift 
D^^'>{x,t) = _D'^'(x) can be transformed via a similarity transformation to a Schrodinger-like equation with time- 
independent potential [J, [3, [E] . Hence the FP equation can be exactly solved if the associated Schrodinger equation 
can be exactly solved. It is also due to this connection that time-independent WKB method can be applied to obtain 
approximate solutions of the FP equation Q . 

When the drift coefficient of the FP equation is time-dependent, the FP equation can still be transformed into 
a time-dependent Schrodinger-like equation, only now the potential in its associated Schrodinger equation is time- 
dependent. We show their connection below. 
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We first define tpix^t) = e'^(^'*^/^^W^(a;, t). Substituting this into the FP equation, we find that ■0 satisfies the 
Schrodinger-like equation: 

where the prime and dot denote the derivatives with respect to x and t, respectively. Eq. Q is the time-dependent 
Schrodinger-like equation associated with the FP equation. Owing to the U term, the potential in the Schrodinger 
equation is time-dependent. 

If ^ can be exactly solved, then the original FP equation is also exactly solved. However, solving the Schrodinger 
equations with time-dependent potentials is generally difficult. In most cases one has to resort to approximate or 
numerical methods. 

In this paper, we would like to show that a direct perturbative approach can be useful to solve the FP equation 
containing a small parameter in the drift potential. 

We shall consider drift potential of the form U{x,t) = X^^^o -^"^"(^' where |A| << 1 is a small parameter. We 
introduce the ansatz ^ — exp{S{x,t, X)/D), where S is expanded in powers of A: 

oo 

5(a:,i,A) = ^A"5„(:r,t). (3) 

n=0 

From ^ we find that S{x,t, /) satisfies 

S = DS" + S'^ + U, (4) 
U ^ ^U"-\U'^ + \U. (5) 

Substituting ([3]) into @ and collecting terms of the same order in A, we arrive at a set of differential equations 
determining the functions Sn- 



So = DS'' + S'^ + ^U'' - \u[? + luo, 



Si = DS'^ + 2Sl^S[ + terms in U of the order A, (7) 
S2 = DS2 + 2S0S2 + 5f + terms in U of the order A^, (8) 



S'„ DSll + S'kS'^_^ + terms in U of the order A", n > 0. (9) 

We note that the probability density Wq ~ cxp{So / D) with Sq satisfying ^ is the solution of the FP equation with 
constant D and drift potential Uo{x,t). 

By solving Sq, Si, S2, ■ ■ ■ from the above equations, we can obtain an approximate solution of ■0, and hence of 
W{x,t). We shall illustrate this by two simple examples below. 



III. FP EQUATIONS WITH DRIFT POTENTIAL U{x,t) = \xV{t) 

Let us consider a class of FP equation with a weak drift potential of the form U{x,t) = XxV{t), where V{t) is some 
finite function of time and |A| << 1. In this case, Uq — Un — {n > 2). This FP equation can thus be viewed as a 
diffusion equation (or Wiener process) perturbed by a weak time-dependent drift potential. 

Let us assume the initial profile of W to be the delta-function, i.e. W{x, t) 5{x) at t = as A — > 0. Then the 
solution of ^ is 

So{x,t) = -^\n{ATTDt)-'^^, (10) 



giving the probability density 



Wo = e o = , exp . (11) 

V ^Dt) ^ ' 
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We have Wodx = 1. This is the well known solution of the diffusion equation having the delta-function as its 
initial profile. 

5*1(2;, i) is determined from 5*0 by ([7]) 



It is easily checked that (IT2|) is solved by 



S,^DS';-^S[ + ^V. (12) 

S^{x,t) = |y(t)-|F(t)+co (13) 

V{t) = J V{t)dt + ci. (14) 

Here cq and ci are arbitrary real constants, cq can be absorbed into the normalization constant of W{x,t). We 
demand that ci be so chosen as to render Si finite in the limit t ^ 0. For instance, if V{t) — cos(ijji), we shall choose 
ci = 0, whilst for V{t) = sin(wt) we have ci = l/oj. 
For S2{x,t), it is determined by ^ 

S, = DS'i + \ {v{t) \v{t)^ "_ml^ (15) 

Solution of (Uni is 



S2ix,t)^-^Vitf+c,^+C2. (16) 

Again, C2 can be absorbed into the normalization constant. We demand that 52 be finite for all a: as t — > 0, hence we 
set C3 — 0. With these forms of ^o, Si and 5*2, one finds from ^ that S3 [n > 3) satisfies the equation: 

S3 = DS'^^^S'„ (17) 

which is solved by S3 — c^x/t + C4 with arbitrary constants C4 and C5. Again finiteness of 6*3 as t — > rules out the 
x/t term, and the constant term can be absorbed into the normalization constant, and hence we set 04 = 05 = 0. 
Continuing this process, we find that all subsequent S'„'s also satisfy (1171) . and hence all Sn must be taken as constant, 
which can all be set to zero. 

Hence, in this example we can solve all the Sn from the set of equations (l7|), ^ to ([9]). The probability density 
W{x^t) is given by the following expression: 

W{x,t) = e-i%V(t)^(So + XS,+X-S,)/D 

^ ^«^pf-77^(^ + ^"^W)'V (18) 



W in (fTS]) is normalized. One notes that is very similar to the solution (|lip of the diffusion equation, only with 
the X in Wq being replaced by a; = x + XV (t). It turns out that W in (fT8| is indeed the exact solution of the FP 
equation with U{x,t) = XxV{t) for any value of A. In fact, in terms of the new variable x and t, the FP equation 
reduces to the diffusion equation: 

dW{x,t) _ ^ d^W{x,t) 

-^r- dx^ ■ ^^^^ 

Eq. (|18p is the exact solution of this diffusion equation with initial profile 6{x + W{Q)). The peak of W{x, t) shifts 
according to the function XV {t)). For instance, if V{x) ~ sin(wt) or cos(a;i), the peak oscillates about a; = as 
W{x,t) spreads outward. 

IV. UHLENBECK-ORNSTEIN PROCESS WITH A SMALL DRIFT POTENTIAL 

Next we consider the Uhlenbeck-Ornstein process [H, Q, which is described by the FP equation with constant 
diffusion coefficient and a time-independent drift potential U{x) = Xx"^ /2. This process is exactly solvable, with the 
probability density given by 



^) = V 2.i^(l-e-2>^) {-m^^^ I • (2°) 
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We would like to apply the above perturbative approach to this process when A is small. Hence we are treating the 
Uhlenbeck-Ornstein process as a perturbed Wiener process. 

The associated Schrodinger equation ([2]) and C/ in ([5]) do not contain U in this case. 5*0 is still to be solved from 
If we assume an initial profile W{x, 0) S{x) as A ^ 0, then 5*0 is again given by PH)) . 

With this 5*0, and §1 read 

S, = DS'{-^S{ + ^, (21) 
^2 = Z?5^'-|5^ + 5f (22) 

Eq. dni) is solved by 

1 X 

Si ^ -Dt + ci— + cq, Co, ci: arbitrary constants. (23) 

As before, finiteness of Si requires ci = 0, and cq is absorbed into normalization constant. Plugging 5*1 = Dt/2 into 
(f22|) . we obtain the solution for 5*2: 

52 = ~Y2^^^ ^ 12^^^ ^ '^''T ^ arbitrary constants. (24) 

By the same reasons given before, we set cz — and leave out C2. Putting all the solutions together, we have, up to 
the A^ terms, the approximate expression for W 

- Tib''"- (-3^ (' + ^' + 5^''') + 1" - T^'"'') - P^' 

Now using the series expansion of ln(l + x) « x — x'^ /2 + . . . for small x, we can rewrite the last two terms in the 
exponent of ([25|) as 

lxt-^Xh'^^\n(^l + Xt + ^Xh^y (26) 

Then (l25l) becomes 



. X /l + At + iA2i2 / a;2 / 1,9,, 



given by ((27)) is normalized. It is easily checked that ([271) is the same as the approximate expression obtained from 
(|20p when the factors in the square-root and the exponent are expanded up to A^ terms. Hence, the result obtained 
by the perturbative approach is consistent with the exact result. 



V. SUMMARY 



In summary, we have applied a direct perturbative theory to solve certain FP equations, which have constant 
diffusion coefficients. Two examples are used to illustrate the method. In the first example the drift coefficient 
depends only on time but not on space. In the second example we treat the Uhlenbeck-Ornstein process with a small 
drift coefficient. These examples demonstrate that such perturbative approach is feasible for obtaining approximate 
solutions of FP equations. 
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